Multi-component gap solitons in spinor Bose-Einstein condensates 
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We model the nonlinear behaviour of spin-1 Bose-Einstein condensates (BECs) with repulsive 
spin-independent interactions and either ferromagnetic or anti-ferromagnetic (polar) spin-dependent 
interactions, loaded into a one-dimensional optical lattice potential. We show that both types of 
BECs exhibit dynamical instabilities and may form spatially localized multi-component structures. 
The localized states of the spinor matter waves take the form of vector gap solitons and self-trapped 
waves that exist only within gaps of the linear Bloch-wave band-gap spectrum. Of special interest 
are the nonlinear localized states that do not exhibit a common spatial density profile shared by all 
condensate components, and consequently cannot be described by the single mode approximation 
(SMA) frequently employed within the framework of the mean-field treatment. We show that the 
non-SMA states can exhibits Josephson-like internal oscillations and self-magnetisation, i.e. intrinsic 
precession of the local spin. Finally, we demonstrate that non-stationary states of a spinor BEC 
in a lattice exhibit coherent undamped spin-mixing dynamics, and that their controlled conversion 
into a stationary state can be achieved by the application of an external magnetic field. 

PACS numbers: 03.75.Lm,03.75.Mn 
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I. INTRODUCTION 



The development of optical trapping techniques for 
Bose-Einstein condensates (BECs) has enabled the con- 
finement of atoms independently of their spin orientation 
in so-called spinor condensates [HQ- Diverse experimen- 
tal and theoretical activities had aimed to characterise 
the effect of the spin degree of freedom on the macro- 
scopic quantum state. A multitude of exciting phenom- 
ena which are not present in a single component BEC 
have been revealed. These include spin textures [1, 0|, 
spin domain formation 0, Q and multi-component soli- 
tons @, H, 0. One of the most frequently used op- 
tical traps for ultra-cold matter is an optical lattice, 
i.e. a standing light wave far detuned from the atomic 
resonance. However, up to now, the effect of the spin de- 
gree of freedom on the key properties of BECs in periodic 
potentials has not been fully explored. 

A unique property arising due to the interplay of the 
inherent nonlinearity of matter waves and the period- 
icity of the optical lattice is the spatial localization of 
a BEC with repulsive inter-atomic interactions. Bragg 
scattering of matter waves with repulsive atomic inter- 
actions can reverse the effective diffraction properties 
of the BEC wavepackets and lead to formation of lo- 
calized states with energy eigenvalues within the gaps 
of the linear Bloch-wave spectrum of the lattice [see 
Fig. QJa,b)]. Two types of such localized structures have 
been studied in the recent years - gap solitons and self- 
trapped states. Single component gap solitons in a one- 
dimensional (ID) lattice have been extensively studied 
theoretically [nl El El El EI El and observed ex- 
perimentally 16 1 . Two-component gap solitons in lat- 
tices have also been investigated [H E3- These non- 



numbers. In contrast, the experimentally observed self- 
trapped states 11911 can confine large number of conden- 
sate atoms fltl |20L l2l| . and can be interpreted as gap 
waves - truncated Bloch states that link gap solitons with 
spatially extended nonlinear Bloch waves in the lattice 

In this paper we present a detailed analysis of the non- 
linear localization of spinor (three-component) BECs in 
the presence of a ID lattice, within the framework of the 
mean-field theory. We consider BEC clouds with either a 
ferromagnetic ( 87 Rb) or polar ( 23 Na) spin-dependent in- 
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FIG. 1: Band-gap structure of the linear Bloch-wave spectrum 
for a periodic potential of a one-dimensional optical lattice 
Vol = Vo sin 2 (x).(a) Variation of bands (shaded) and gaps 
(clear areas) with the potential depth, fJ.(Vo), and (b) eigen- 
value spectrum fi(k) reduced to the first Brillouin zone shown 
for Vo = 2. Points (c) and (d) correspond to the ground and 
excited state linear Bloch waves with the typical spatial den- 
sity profiles in a single hyperfine component of the condensate 
shown in panels (c) and (d), respectively. 
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teraction, and show that the lattice enables the existence 
of multi-component localized states of both types of the 
spinor condensates. We show that the spinor states in 
both of the two condensate types may be localized as 
vector gap solitons that can exhibit cither ferromagnetic 
or polar spin structures regardless of the type of spin- 
dependent interaction. In addition, we find and charac- 
terise the multi-component self-trapped gap states that 
naturally extend the concept of single-component self- 
trapped gap waves [20j | . 

Both the spinor gap solitons and gap waves can be 
adequately described within the single-mode approxima- 
tion (SMA) which assumes that all components of the 
composite localized state share the same spatial profile. 
However, we also reveal that both gap solitons and gap 
wave states may exist simultaneously in different spin 
components forming a hybrid localized state, which no 
longer satisfies the SMA. Such composite states may ex- 
hibit effects of self-magnetisation and internal Josephson- 
like dynamics [Hj]. 

Furthermore, we consider non-stationary dynamics of 
the spinor gap solitons and examine two distinct meth- 
ods to induce coherent spin-mixing oscillations, through a 
change in the relative phase of the condensate, or through 
a population imbalance in a stationary state configura- 
tion. We find that the presence of the lattice potential 
prevents fragmentation of a spinor BEC and damping of 
the coherent spin-mixing evolution. Finally, we demon- 
strate that application of a magnetic field [231 ] can arrest 
the spin mixing, which lays the foundations for controlled 
and flexible state manipulation. 

The paper is organized as follows: section [III where we 
introduce our model is followed by section IIIII in which 
basic properties of stationary states are analysed. In sec- 
tion IIVI we examine the dynamical instability of spinor 
Bloch waves and spinor matter wave localization. Sec- 
tions |Vl I VII and IVIII concern distinct types of localized 
stationary states. And finally, in section IVIIII coher- 
ent spin-mixing dynamics in multi-component localized 
states is discussed. 



II. MODEL 

In the framework of mean-field theory, the wavefunc- 
tions, = v f( r )±i,o, of the three hypcrfinc com- 

ponents (mp = ±1,0) of an F = 1 spinor condensate 
obey the following set of coupled Gross-Pitaevskii equa- 
tions Q: 

i^i = C9±i + X a (|*±i| 2 + |*o| 2 - |* T i| 2 ) *±1 

i^ = £*o + A a (|tf_i| 2 + |<M 2 )*o (1) 
-r-2A„* + i*_i*o, 





a [a B \ 


a2 [asj 


A s 


A a 


8/ Rb 

23 Na 


101.8 
50.0 


100.4 
55.0 


1.49 x 10 _2 
4.08 x 10 -3 


-6.94 x 10~ b 
1.28 x 10~ 4 



TABLE I: Scattering lengths a , a 2 for 87 Rb and 
23 Na [29l ] atoms given in units of the Bohr radius as- The 
dimensionless ID coupling constants A s , A a used throughout 
the paper are obtained according to Eqs. ([2} (see text). 

where L = - ±V 2 + V(r) + X s (|*_i| 2 + |* | 2 + |*+i| 2 ). 
We assume the confining potential potential to be of 
the form V(r) = ^J\(y 2 + z 2 ) + V l{x), where ujj_ is 
the trapping frequency in the direction perpendicular 
to the lattice, and Vol = Vo sin 2 (7rx/c?) is a periodic 
potential created by a ID optical lattice, with the pe- 
riod d. The wavefunctions are normalised to the num- 
ber of atoms: J dr \^j\ 2 = Nj, where j = ±1,0 and 
SiVj = J\f is the total number of atoms in a conden- 
sate. Equations U) were made dimensionless by measur- 
ing length, time and energy in units: aj_ = (H/mu!±)^ , 
wj 1 , and fku±, respectively. In this work we consider 
both condensates of 87 Rb and 23 Na atoms, and assume 
strongly anisotropic clouds with transverse trapping fre- 
quency uj± = 2tt x 230 Hz 3> uiu (25[. The transverse 
oscillator length is then a± = 0.71 /im for 87 Rb atoms 
and a± = 1.38 /xm for 23 Na atoms. In what follows, we 
set d = 7r, which for our parameters corresponds to the 
lattice spacing of 2.23 fim for 87 Rb and 4.34 /jm for 23 Na 
atoms. 

Spinor condensates exhibit both spin-independent and 
spin-dependent nonlinear inter-atomic interactions, with 
the sign of the spin-dependent mean-field energy effec- 
tively determining the properties of the ground and ex- 
cited states. A BEC of bosonic spin-1 atoms can be ei- 
ther ferromagnetic (e.g. 87 Rb) or polar (e.g. 23 Na) in na- 
ture [26l . [2?| , with the two types of BECs having different 
average values of the total spin and different symmetry 
properties. The coupling coefficients for "symmetric" 
spin-independent and "antisymmetric" spin-dependent 
interaction terms in Eqs.{T]) are: 

Here ao and 02 characterise the s-wave scattering of two 
atoms in the channels with total spin equal to and 2 re- 
spectively. The measured values of the scattering lengths 
for 87 Rb (see tablcQJ result in A Q < [28[ and a. ferromag- 
netic ground state of the condensate. In contrast, a 23 Na 
BEC has A a > and the condensate is polar [29|. The 
strong inequality | A a | <C | X s | , which holds for both types 
of the BECs, means that the spin-dependent interactions 
are much weaker than the spin-independent ones. 

Under the assumption of a strong anisotropy and a 
tight confinement in the dimensions transverse to the op- 
tical lattice, the model ((T|) can be further reduced to a 
system of mean-field equations for the axial components 
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of the wavefunctions. By introducing separable wave- 
functions, = ipj(x)')p±(y, z), with the transverse com- 
ponents ip±(y,z) determined by the ground state of the 
tight harmonic trap (30j . and integrating Eqs. ([T]) over 
the transverse directions {y, z}, the model for the spinor 
condensate is reduced to the ID evolution equations for 
%j)j(x). The coupling constants arc renormalised to values 
listed in table [J For later reference, we define the total 
density of atoms as nxot = £ rij(x,t) = £ \tpj{x,t)\ 2 , 
where n±i.o are the atomic densities of different hyper- 
fine states. 

The reduced dynamical system of ID GP equations for 
is characterised by the conserved total number of 

atoms J\f = J^^dx (n_i + n + \ + no) and Hamiltonian 

(total energy of the system), H = f_°° h(x)dx, with the 
Hamiltonian density: 



h = E ( \ ^ 



Vol nj + — n 



(3) 



J=0,±1 
j,k=0,±l 

The spin-dependent part of the Hamiltonian density, 
h a = (A a /2) |sj 2 , is proportional to the square of the 
local spin expectation vector: 



s = V-' T FV.' = n T otC T FC 



(4) 



Here ip = (VM-Ij V'Oj V ; -i) T is a three-component wave 
function, £ = (C+i, Co 5 C-i) T i s the three component 
spinor, normalised according to the condition C* T C = 1; 
and F = (F x ,F y ,F z ) T is the vector of spin opera- 
tors [1, HH, so that the local spin expectation vector 
yields: 



y/2 Tle{^4'+i + ip-iM 



(5) 



n_i 



In what follows, we will only refer to the normalized spin 
expectation vector (local spin) 



S = s/n To t, 



(6) 



where we note S = (S x ,S y ,S z ) T . Furthermore, we define 
the averaged spin expectation vector (i.e. total spin) as: 



(S) = / dx s(x)/ / dx \i/j(x)\ 



(7) 



In the definition ([7]) the average spin is normalized such 
that (S) takes the maximal value 1. In the absence of a 
magnetic field, due to the rotational symmetry the total 
angular momentum of the system (S) is conserved. We 
note that (S z ) = f_°° {n~i — n + i) dx, known as the sys- 
tem's magnetisation, is the difference in the fraction of 



atoms in the Zeeman states mp = ±1. The conservation 
of the magnetisation in this system has a deep ana logy 
in the optical parametric four- wave mixing process [3 11 ] - 
The four-wave mixing interaction, involving e.g. the anni- 
hilation of a pair of atoms with mp = +1 and mp = — 1, 
and creation of two atoms in the mp = state, leads to 
phase conjugation in quantum optics and to matter-wave 
phase conjugation in BECs {32} . 



III. GENERAL PROPERTIES OF STATIONARY 
STATES 

The stationary states of a spinor condensate can be 
found in the form ipj(x,t) = 4>j(x) exp(— i[ijt) exp(i9j), 
where </>j(a;) are real wavefunctions, \ij are chemical po- 
tentials, and 8j are absolute phases of the three hypcrfine 
components. As usual for physical systems with paramet- 
ric interactions 31| , the phase matching condition for the 
chemical potentials, 2/io = (J>+i + must be satisfied 
for a stationary state to exist. In addition, the relative 
phase between the components, A6* = 28q — (8 + i + 6Li), 
can only take values or ir. The spatial profiles of the 
stationary state wavefunctions, <j)j(x), are then set by the 
following coupled time-independent ID equations: 

fJ—i<j)-i = £id0-i + A a ((j) 2 _ 1 + 4>l — 4> 2 +1 ) </>_i 



HO<Po = £lD0O + A a (<^_j + </> +1 ) (j) (8) 

± 2A a + i0_i0o 
M+i0+i = Ad0+i + A a {4> 2 +1 + (j>l - (fri)) + i 
± A Q 0o0_i 



where C 1D = -\-§^ + Vol (a) + A s +$ + <t> 2 +1 ). 
The plus sign in front of the last term in Eqs. ([5]) corre- 
sponds to stationary states with all three components in- 
phase with one another and therefore A9 = 0. The minus 
sign is realized for A8 = ir, and then one of the compo- 
nents (mp = +1 or mp = —1) of a stationary state is 
out-of-phase with the others. We note that both types of 
stationary states always exist in the system, regardless of 
the sign of A a . However, it is easy to see @, that the in- 
phase stationary states minimise the spin-dependent part 
of the total energy ([3]) of the system, H a = f^°° h a [x)dx, 
for A a < 0, while the out-of-phase stationary states min- 
imise the spin-dependent energy for A a > [26]. The 
spin-dependent part of the Hamiltonian thus determines 
the ground state of the condensate: ferromagnetic for 
A Q < 0, and polar for A a > @, H, [H, SI . 

For large atomic densities and atom numbers the spin- 
dependent interaction energy H a is sufficient to create 
spatial spin textures in lattice-free condensates and thus 
leads to peculiar miscibility properties of spinor BECs [f| • 
Due to the strong inequality of the coupling constants 
|A a | <C | A s |, the spin-dependent interaction has little ef- 
fect on the spatial density distributions of the station- 
ary states for low and moderate peak densities. Hence 
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in small atomic ensembles the internal spatial dynamics 
is suppressed [H, HH HH , although H a still determines 
the phase structure of the ground state and governs spin- 
mixing dynamics [34| . The spatial density structure of 
the condensate can be determined from Eqs. |(5J). 

When the size of a condens ate is sm aller than the spin 
healing length £ sp i n = 2ir/ y / 2\X a \n. the three compo- 
nents of the condensate can be treated as if trapped in 
the same effective potential created by an optical lat- 
tice and the combined mean-field densities of the spinor 
components. This single mode approximation has been 
widely employed to describe the properties of the ground 
eigenstates of this effective potential [HI, [H, [H, l37| . 
Within this approximation, a stationary vector wavefunc- 
tion of a spinor condensate can be represented in the 
form: 4>j(x) = <p(x)Q, and all the components share the 
common density profile cf> 2 that obeys the single Gross- 
Pitaevskii equation: 



ld 2 4> 
'^dx 2 



^ol(x)0 + A s 3 , 



(9) 



where fi = fiQ = = The spatial distribution of 

the local spin S is governed by the relative distributions 
of condensate components. As a consequence of a shared 
density, the distribution of the local spin is reduced to a 
single point in the 3D spin space. 

Within the SMA, the in-phase (ferromagnetic) station- 
ary states are characterised by the spinor [2(| : 




cos 2 (/3/2) 
/ 2cos(/3/2) sin(/3/2) 
sin 2 (/3/2) 



(10) 







t~fr sin/3 




-I 


cos/3 






v 75 sin/3 



where (3 is an arbitrary rotation angle around the z-axis 
(see Fig. 2] for the definition). Accordingly, the squared 
expectation value of spin in this state is (S) 2 = 1. The 
spinor for out-of-phase (polar) stationary states is given 
by M- 



(11) 



and (S) 2 = for this state. Our numerical solutions of 
Eqs. ©, which are presented in sections HVl M and IVI1 
confirm the validity of the SMA for a variety of ex- 
tended and spatially localized nonlinear stationary states 
(of both polar and ferromagnetic type) for spinor BECs 
in optical lattices. However, as will be shown below, there 
exist nonlinear stationary states that cannot be described 
by the SMA. 



IV. STABILITY OF SPINOR BLOCH WAVES 

When the interatomic interaction can be neglected the 
eigenstates of Eqs. © take the form of linear Bloch 



states: cj>(x) ~ (f>k(x)exp(ikx) , where <pk{x) = <fik{x + d) 
has the periodicity of the lattice, and the wave num- 
ber k is selected from the first Brillouin zone (BZ), 
i.e. —1 < k < 1 in our units. Then the linear spec- 
trum /j(k) of the Bloch states has a characteristic band- 
gap structure with the lowest two bands shown in Fig. [T] 
In the presence of interactions, the most fundamental 
stationary states of a BEC in an optical lattice are the 
periodic, spatially extended nonlinear Bloch waves [ll[ 
that transform into linear Bloch waves for sufficiently low 
atomic densities (nonlincarity). 

It is known that Bloch states of a single-component 
BEC with repulsive atomic interactions, corresponding 
to the top edge of the ground band (i.e. the edge of the 
first Brillouin zone), exhibit complex nonlinear dynamics 
due to modulational instability (MI) [3o |. This instabil- 
ity arises because the anomalous diffraction at the band 
edge effectively changes the sign of the atomic interac- 
tions and can lead to effects that are normally absent 
in a repulsive condensate, such as formation of spatially 
localized (bright) gap solitons. As a result of such in- 
stability development, gap solitons trains can form in a 
process of condensate evolution at the edge of a Brilllouin 
zone SSllii. 
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FIG. 2: Development of modulation instability for F = 1 
spinor Bloch states in (a) ferromagnetic ( 87 Rb) and (c) polar 
( 23 Na) BECs at the top edge of the first band (Vo = 4, p = 
1.4). Shown are the peak densities of the tuf — +1 (dashed 
line), mF ~ (dash-dotted line) and m,F = — 1 (dotted line) 
condensate components at x = 100. For comparison, the 
behaviour of the modulationally stable ground state nonlinear 
Bloch waves at /i = 1.28 are shown for a (b) ferromagnetic 
( 87 Rb) and (d) polar ( 23 Na) BECs. 

In the absense of a lattice, the stability properties of a 
spinor BEC confined in a harmonic trap are well estab- 
lished. Namely, a condensate with ferromagnetic type of 
spin-dependent atomic interactions is found to be modu- 
lationally unstable, and the polar BEC is stable [37l.[4l|. 
Here, by performing nymerical analysis of condensate 
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FIG. 3: Instability-induced dynamics for a 23 Na condensate 
at the top edge of the ground band (fi = 1.12, Vo = 2). Shown 
are snapshots of the total BEC density (solid lines) with the 
periodic OL potential (gray shading, magnified by a facto of 
20) at the times (a) t = 0, (b) t = 20, (c) t = 40, (d) t = 100, 
(e) t = 120, and (f) t = 200. 



evolution [521 ]. we find that Bloch states corresponding 
to the upper edge of the ground band [see Fig. 1(c)] 
for both, ferromagnetic ( 87 Rb) and polar ( 23 Na) BECs 
are modulationally unstable. In contrast, ground Bloch 
states corresponding the lower edge of the band [see Fig. 
1(d)] for both 87 Rb and 23 Na are modulationally stable, 
irrespective of the phase relationahsip between the hy- 
perfme components. 

Typical results of our numerical simulations demon- 
strating the stable and unstable dynamics of Bloch states 
in an optical lattice of moderate depth are presented 
in Fig. [21 For the chosen parameters, the duration of 
evolution sufficient to detect the onset of MI is 200 time 
units (140ms). The onset of the MI occurs at t ~ 130 
(90 ms) for the unstable states of both ferromagnetic and 
polar condensate types at the top of the first band, as 
seen in Fig. [5Ja,c). The time scale of the MI develop- 
ment is comparable with the time scale associated with 
density-dependent interactions, t s = A s riTot- The peak 
densities of the initial states nxot = 8.9 for 87 Rb and 
n Tot = 32.7 for 23 Na give t s = 0.134 (92.5 (is) for rubid- 
ium and t s = 0.133 (92 /is) for sodium. For the same 
lattice parameters, Bloch states from the bottom of the 
BZ demonstrate stable dynamics, as shown in Fig.[5(b,d). 

The dynamical instability of the Bloch states at the 
band edge underpins the existence of the spatially local- 
ized states. In order to further investigate the dynam- 
ics of the Bloch states beyond the onset of the MI and 
initiate formation of localized structures [40j . we intro- 
duce initial periodic phase modulation, cxp{i cos(<5x)}, 
imprinted onto all the hyperfinc components in the ini- 
tial state. As an initial state, we take an unstable Bloch 
wave corresponding to the top of the first band in a rela- 
tively shallow (Vo = 2) lattice. This results in larger peak 




FIG. 4: Examples of spatial structure of the (a,c,d) 
ferromagnetic-type (in-phase) and (b) polar-type, (out-of- 
phase) gap solitons of a polar 23 Na condensate in an optical 
lattice (Vo — 4, fi = 2.8). Red, green, and blue lines show the 
4>-i, 4>o, and 4>+i components, respectively. The correspond- 
ing values of the local spin vector S = (S x ,S y ,S z ) are plotted 
in the plane S y = of the 3D spin sphere (shaded); /3 is the 
Euler rotation angle. 



densities nxot of 16.6 for 87 Rb and 61.7 for 23 Na, and 
shorter interaction time scales t s ~ 0.25 (0.2 ms) for both 
atomic species. In Fig. [3] we show snapshots of the total 
density riTot in a shallow lattice for anti-fcrromagnctic 
matter waves. The wavelength of the phase modulation, 
A = 2tv/6, corresponds to ten lattice sites for 5 = 0.2. 
As a result of the phase imprinting, strong density mod- 
ulations develop at the initial stages of the evolution. 
Periodic density modulations trigger modulational insta- 
bility and greatly accelerate the process of the conden- 
sate localization, as demonstrated previously for a single- 
component BEC |40j . 

The scenario of dynamical localization is almost identi- 
cal for both 87 Rb and 23 Na condensates (shown in Fig. [3] 
for a sodium BEC). This dynamics is primarily deter- 
mined by the spin-independent interactions and hence is 
similar to the dynamics of a single component conden- 
sate with repulsive inter-atomic interactions [3^ ]. The 
development of the localized structures is complex, with a 
transition from single-peak soliton-like structures [Fig. [3] 
(b)] to unstable "bound states" of two neighbouring soli- 
tons [Fig. [21(d)]. The nature of the spin-dependent inter- 
actions turns out to play a role in the interaction of local- 
ized structures. Namely, for polar ( 23 Na) condensate the 
neighbouring localized states tend to repel each other due 
to the presence of out-of-phase components. Conversely, 
for a ferromagnetic ( 87 Rb) condensate, the neighbour- 
ing solitons with all three components in-phase attract 
each other. In general the "bound states" do not persist 
over long evolution times and the dynamically formed 
solitonic structures eventually disappear [Fig. [3(f)]. 



V. SPINOR GAP SOLITONS 

Unlike the spatially extended nonlinear Bloch waves, 
spatially localized states of a condensate in an optical 
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lattice can exist only within the spectral band gap. Here 
we describe the properties of stationary three- component 
(vector) gap solitons. To this end, we numerically solve 
time- independent equations i[8]). for both types of spin- 
dependent nonlinear interactions, looking for spatially lo- 
calized solutions within the gap. The results of our find- 
ings are summarised in Fig. [4j where we plot the local 
spin vector in the {S X ,S Z } coordinate plane (S y = due 
to the phase matching condition for stationary states). 
As shown in the figure, all stationary three-component 
gap solitons can be divided into the two types, in- 
phase (ferromagnetic-type) and out-of-phase (polar-type) 
states. We stress that both types of solution exist for 
both ferromagnetic (A a < 0) and polar (A a > 0) types 
of spin-dependent interactions. In agreement with the 
single mode approximation, all the in-phase states can 
be described as a point on a circle of unit radius on 
the {S X ,S Z } coordinate plane. The peak density of the 
mp = component ranges from at point (d) , to half of 
the total peak density, n-Tot/2, at point (c). In contrast, 
all the out-of-phase gap states correspond to the point at 
the origin {S x , S z } — {0, 0}. 

Profiles corresponding to point (a) in Fig. [5] represent 
a typical three-component in-phase gap soliton, where 
the strong localization occurring in a nonzero spin com- 
ponent, e.g. m p = —1, induces the localization in the 
other two hyperfine states. In contrast, for the in-phase 
gap solitons corresponding to S z = and represented by 
point (c), the amplitudes of the component wavefunctions 
for mp = +1 and mp = —1 are equal, see Fig. Etc), so 
that the populations of the states are: no = nxot/2 and 
n_i = n + \ = rtTot/4. Two of the spinor components 
are equal: £+i = C-i- At point (d) (or its diametrically 
opposite), Co = and two spin components have zero 
amplitudes to yield a maximum (or minimum) value of 
the magnetisation. The example in Fig. [UJd) shows a 
gap soliton with all atoms in the mp = —1 state, which 
maximises the magnetisation. 

All the out-of-phase solutions to Eqs. © have 4>+\ and 
<^>_i components that are equal in their absolute values 
but 7r out of phase, which effectively leads to a cancella- 
tion of the spin interaction terms in model |(5J) and a zero 
value of the magnetisation. It should be noted that some 
states represented by points (b) and the in-phase state 
indicated by point (c) in Fig. 0] can not be distinguished 
in an experiment by simple reconstruction of the com- 
ponent density profiles, as they might differ from each 
other only by a relative phase. The recently developed 
nondestructive method of polarization-dependent phase- 
contrast imaging 42j opens up new possibilities for prob- 
ing spinor BECs, and can ultimately enable the practi- 
cal identification of the in-phase and out-of-phase states 
and lead to recognition of ground states and long-lived 
metastable spin states. 

In general, stationary gap solitons of the in-phase type 
can have atoms populating one or three hyperfine states, 
while the solitons of the out-of-phase type can also form 
two component states. In Fig. [5] we show the variation of 




FIG. 5: Comparison of the values of gap soliton peak ampli- 
tudes, (j>+i and (j>-i at x — (dashed lines), with the ampli- 
tudes of a homogeneous condensate (solid lines) for (a) out-of- 
phase and (b) in-phase stationary states of a S7 Rb condensate 
in a lattice with Vq = 4, plotted as a function of cj>Q. Gap soli- 
ton is taken at n = 2.8, and the homogeneous condensate at 
/i e fi = 2.0, to account for the shift of the spectrum from zero 
(see text). Analytically obtained parameters of the semi-axes 
of the ellipses and the ellipse displacement are (a) A = 11.59, 
B = 8.19 and (b) A = 8.21, B = 5.8 and A+i = B. 



the soliton peak amplitudes in different hyperfine states, 
depending on a given peak amplitude in the mp = 
state, for a particular set of lattice parameters and a 
chemical potential within the first gap. We compare 
these with the peak amplitudes calculated for a homo- 
geneous condensate. The introduction of a periodic po- 
tential shifts the chemical potential of the extended so- 
lution with the lowest chemical potential (the bottom of 
the first band) relative to the homogeneous condensate 
(the solution at Vq = 0), as can be seen in Fig. [TJa,b). 
Thus in comparing peak amplitudes it is important that 
the same effective chemical potential is used. To this end 
we use for the homogeneous condensate ^ e ff = A* — f-o 
where is the chemical potential of the bottom of the 
first band. 

Considering the out-of-phase states in a homogeneous 
case we set, / = = — V^-i and g — <f>o, so 

that Eqs. ijHJ) simplify to: 

Meff / = K (/ 2 + 9 2 ) /; Mcff .9 = A Q (f 2 + g 2 ) g. 

(12) 

This set of equations has an implicit analytical solution 
which can be written in the form of an ellipse equation 
for <A»: 



A 2 B 2 



1. 



(13) 



where A = y' /j, e g / X s and B = /i e ff/(2A s ) are the semi- 
major and the semi-minor axes of the ellipse respectively. 
For the amplitudes of an in-phase state, equations 
cannot be simplified, and give the implicit solution: 



0§ 

.4 2 



A+i) 2 



B 2 



1. 



(14) 



Here the ellipse parameters A, B, and the dis- 
placement of the ellipse centre A+i are defined as: 
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FIG. 6: Self-trapped gap states at fi = 2.8 corresponding 
to the nonlinear Bloch waves of a 87 Rb condensate (a) in the 
ground state and (b) at the edge of the first band of the optical 
lattice potential (shaded) with Vo = 4. Red, green, and blue 
lines show <j)-i, <j)o an d cj>+i hyperfine states, respectively. 



A 2 = [2(\ s + X a )/^ eS }-\ B 2 = [4(A S + X a )/fx eS ] 1 and 
A+i = h\/ /i c ff/(^s + A a ). These solutions, derived from 
the lattice-free equations, predict the amplitudes of the 
multi-component gap solitons with good accuracy (see 
Fig. © . 

As both the in-phase and out-of-phase types of states 
exist in any spin-1 BEC, the generation of both classes 
of gap solitons could be experimentally feasible. A small 
atom number BEC cloud of optically trapped atoms can 
initially be transferred to the desired mj? = +lormp = 
— 1 spin state by Landau-Zener rf-sweeps [43[. Then a 
single-component gap soliton (correspondin g to point (d) 
in Fig. [4]) in an optical lattice can be formed [l6[ . Finally, 
by applying a magnetic field perpendicular to the z-axis, 
a spin-rotation can be induced to change relative popu- 
lation of the spin components within a stationary state. 



FIG. 7: Non-SMA stationary states for a Na condensate: 
(a) centred on the potential maximum (i.e. off-site) and (b) 
centred on the potential minimum (i.e. on-site). Both station- 
ary states are found for the values of the chemical potentials: 
fi+i = 2.8, fi-i = 2.7 and fio = §(/i+i + M-i) = 2.75. Red, 
blue, green lines mark <f>+\, <f>-i and <j>o components, respec- 
tively. Shaed areas show the corresponding optical lattice 
potential (Vo = 4). 



suits in possibilities for different relative populations of 
the hyperfine components, similar to the case of spinor 
gap solitons (as shown in Fig. [4|) . In addition, they also 
form in-phase and out-of-phase spin configurations for 
any type of spin-dependent interactions. Examples of 
out-of-phase truncated nonlinear Bloch waves for 87 Rb 
are shown in Fig. [?5J 

We note, that the spatial structure of self-trapped gap 
waves can be adequately described within the single- 
mode approximation, and hence Eqs. (fTU)l and (fTTj) de- 
termine relative fractional populations for different ori- 
entations of spin in a self-trapped state. 



VI. SELF-TRAPPED STATES 



VII. BEYOND THE SINGLE MODE 
APPROXIMATION 



A second type of spatially localized states existing in a 
lattice, the "self-trapped" states, have recently been in- 
vestigated in a single-component condensate 0, HH • 
These states are nonlinear Bloch waves which at suffi- 
cient densities experience de-tuning of the chemical po- 
tential deep into the linear band gap. Consequently they 
may be localized by Bragg reflection and "truncated" 
spatially. The localization has the same physical origin 
as for the gap solitons, and therefore these states connect 
the two fundamental types of nonlinear lattice states, the 
extended nonlinear Bloch waves and the localized gap 
solitons. The truncated nonlinear Bloch waves may be 
originate from any periodic Bloch state, such as e.g. the 
ground state nonlinear Bloch wave. 

Here we generalize this concept to the multi- 
component Bose-Einstein condensates and demonstrate 
the existence of spinor gap waves (sec Fig. |TJ]). Akin 
to a single component self-trapped state, the ground 
state spinor gap waves only exist beyond a critical lat- 
tice depth, sufficient for localization to take place. Be- 
yond this point the gap waves may be localized with an 
arbitrary spatial extent [2(J. Nonetheless, in contrast to 
a single component BEC the spin degree of freedom re- 



Apart from the "fundamental" gap states considered 
above, and adequately described by the single-mode ap- 
proximation, an optical lattice supports the existence of 
other spinor gap states with more complex spatial struc- 
ture. Under special conditions, off-site and on-site gap 
solitons that are extended over multiple lattice sites and 
centered on a maximum or minimum of the lattice po- 
tential respectively can be found. These and other "ex- 
cited" gap states contain a larger number of atoms than 
the fundamental solitons. Most importantly, due to the 
dramatic differences between the spatial density distri- 
butions in the three hyperfine states, they do not obey 
the single mode approximation. Examples of non-SMA 
states obtained from equations ([5| for a 23 Na condensate 
are shown in Fig. [7J 

The existence of the states that do not obey the SMA 
is counter-inituitive as the size of gap solitons is an or- 
der of magnitude smaller than a typical spin domain [43| . 
These states may be regarded as a composite of the self- 
trapped states and gap solitons, in which a gap soliton is 
localized in the combined effective potential of the self- 
trapped state and the lattice. Straightforwardly, the lo- 
calized structure in nip = is induced by a 7r phase jump 
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2000 



FIG. 8: Spin dynamics of unstable off-site non-SMA state of 
23 Na for parameters given in Fig [9] The magnitude of S x is 
taken at x = §7r. Initially, the local spin oscillates with period 
T = fj- ~ 104.7. For a time interval of nine oscillations, which 
correspond to t = 942.5, the dynamics is stable. Afterwards 
the unstable evolution commences. 



in the density envelope of a strongly populated compo- 
nent, either mp = —1 or mp = +1. Importantly, here 
the spin degeneracy has been removed due to the dif- 
ference in chemical potentials, i.e. fi—i M+i? so the 
populations in each hyperfine component are no longer 
governed by simple rotation through spin space [see Eqs. 
([TO]) and (fTTj) ]. Furthermore, it is clear that within these 
states one component is weakly populated, leading to an 
effective two-component system. Despite this reduction 
in the number of the hyperfine states, the spin becomes 
here an important quantity, with the solutions possessing 
a complex spin structure. 

The spatially separated wavefunctions of the differ- 
ent hyperfine components give V£ ^ and lead to a 
violation of the single mode approximation. Hence the 
spin-dependent part of the total energy H a is no longer 
constant. As an effect the values of the local spin tend 
to cover the whole range < |«S| < 1 and gap solitons 
acquire spin structure across the localization region, as 
shown in Fig. [TU] 

By numerical simulations we determine that the on-site 
non-SMA soliton shown in Fig. [7Jb) is stable. However, 
due to the difference in the chemical potential between 
the components the local spin exhibits spin precession 
dynamics as given by the expression: 



5 x (a;,0)cos (Sfj.t + 66) 
S(x,t)=\ Sy(x,0)sm(Sfj,t + 86) 
S z {x,0) 



(15) 



where 86 = — an d Sfj, = — M+i) ac ~ 

cording to the phase-matching conditions for stationary 
states. For real stationary states the local spin oscillates 
with the frequency w = 8fi. As the non-SMA states are 
made of "twisted mode" gap waves in mj? = ±1 compo- 
nents which are 7r-out of phase from each other as shown 
in Fig. [Jj 86 = it. We have checked that in the regimes 
of their stability the non-SMA states do indeed exhibit 
spin oscillations at the frequency 8[i. We display those 
in Fig. [8] showing S x for the unstable state in Fig. UJa). 
Due to the precession of the local spin in the {S x ,S y } 
plane, each point of the soliton behaves as if it was ex- 
hibiting a permanent magnetic field along the z-axis. 
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FIG. 9: Non-SMA stationary states of 23 Na for parameters: 
jU+i = 1.28, jti-i = 1.4 and /xo = 5(^+1 + = 1-34 and 

Vo = 2. All hyperfine components of the unstable off-site 
soliton at t — (a,b,c) and at t — 2000 (d,e,f) are compared 
against the corresponding profiles of a stable on-site non-SMA 
state (g,h,i). 



Therefore the non-SMA solitons can be thought of as 
being self-magnetised. It is conceivable that polarization- 
dependent phase-contrast imaging [42| could resolve this 
effect. We note that although the local spin exhibits pre- 
cession around the iS 2 axis, the total angular momentum 
is conserved. 



Although the on-site state is stable, the off-site state 
shown in Fig. UJa) is unstable and its stability proper- 
ties depend on the depth of the lattice potential and the 
degree of localization due to the chemical potential. For 
the chosen parameters, this state decays after nine pe- 
riods of local spin oscillation. Due to the self-trapping 
effect, this unstable off-site state does not lose atoms by 
radiation but instead starts to exhibit spin interaction 
induced Josephson-like oscillations [22| between the left- 
hand side and the right-hand side potential wells. To gain 
some insight into the decay mechanism of the unstable 
state, we analyse the development of the instability. Ac- 
cordingly, we plot in Fig. [5] profiles of different mp com- 
ponents of the off-site state at t = (left column) and 
t = 2000 (middle column) corresponding to the initial 
and final time samples of our simulations. The matching 
spin projections at these times are presented in Fig. 1101 
We compare the unstable off-site state against the sta- 
ble on-site non-SMA state (right column in both Fig. [5] 
and Fig. IT0|) . In the course of the unstable evolution, the 
off-site state transforms into an on-site state and then 
evolves into an unstable "bound-state" , similar to those 
emerging from the MI- induced dynamics (see section H*V| . 
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FIG. 10: Projections of local spin of the non-SMA states of 
23 Na (shown in Fig. [9} onto the planes S y = (a,b,c) and 
S x — (d,e,f) of the spin sphere (shaded). Distribution of 
the local spin vector of the unstable off-site non-SMA state in 
its (a,d) initial (a,d) and (b,e) final, i.e. at t = 2000 configura- 
tions. For comparison, panels (c,f) show the spin distribution 
for the stable on-site state. 




FIG. 11: Dynamics of non-equilibrium states in 87 Rb (fi = 
2.8, Vo — 4). (a) Spin mixing oscillations of the populations 
between tuf — —1 (red), rriF = +1 (blue, dashed-dotted) 
and rriF = (black) hyperfine components, (b) The dynam- 
ical trajectories in the no (AS) phase space for the scenario 
presented in panel (a). Initial (final) location in the phase 
space is marked by the green (blue) point. The arrow indi- 
cates the direction of the spin mixing oscillations in the phase 
space. 



VIII. NON-STATIONARY DYNAMICS 

If a spinor BEC is initially prepared in a non-stationary 
spin co nfig uration, a general spin-mixing dynamics oc- 
curs @, [23, [H, HI] , i.e. exchange of population between 
the central and sideband fields. The role of the ferro- 
magnetic (or polar) interactions in this dynamics is to 
raise (lower) the energy of the mp — component rel- 
ative to the energies of the other two components, until 
the system reaches its ferromagnetic or (polar) ground 
state. Predicted theoretically 0, [45j|, the spin-mixing 
dynamics of large-size harmonically trapped condensates 
has recently been experimentally demonstrated @, E(|. 
Unfortunately, the spin-mixing oscillations were observed 
for only a few periods, due to the damping by fragmen- 
tation of the spinor BEC @, Hi]. 

In the framework of the mean-field theory, the co- 
herent evolution of spinor components does not lead to 
damping of spin oscillations. As the size of a single gap 
soliton is much smaller than the spin healing length, in 
contrast to larger condensate clouds 0, |46j], their evo- 
lution can be very accurately described by the single 
mode approximation. The conservative non-equilibrium 
dynamics within the SMA generally consists of peri- 
odic exchange of populations between different hyperfine 
components, which results in the dynamics of spinors: 
(jit) = x /n^)exp(-i^-(f)) [H[. Within the SMA 
the spin-dependent energy landscape is governed by the 
Hamiltonian (36| : 

H = 7Po (l - Po + s/i 1 ~ Po) 2 - m 2 ) cos(A0) + |m 2 

(16) 

with po = no/uTot being the relative population of the 
mp = state. Here m = n +1 — n_i is the value of (the 
local) magnetisation and 7 = A Q Affdx (nxot) 2 is the 
coupling. A8 is the relative phase of the condensate, as 



defined in section IIIII These values arc extracted from 
the multi-component gap soliton. 

We note here that there are two ways to initiate 
the spin-mixing dynamics. The first one is via a non- 
stationary relative phase [HI, E3|. The second is by means 
of population imbalance. Both these procedures give 
equivalent results. It has previously been shown fiil |45|| 
that if the population is averaged over many realizations 
with different values of the initial condensate phase, the 
spin-mixing dynamics can be halted in the averaged pic- 
ture, and the system may reach a steady-state. 

We monitor the relative phases of solitons with the 
introduced population imbalance and plot their trajec- 
tories in the no(A0) phase space, as shown in Fig. QT] 
In the figure we also compare the resulting trajectories 
with the energy landscapes, Eq. (fPo]) . determined sep- 
arately for each scenario. We prepare initial states by 
shifting the relative phases or by introducing a popula- 
tion imbalance in a stationary gap soliton, e.g. as shown 
in Fig. [4] An example of the dynamics in the first case 
is shown in Fig. 111! with the green point in Fig. Illf b) 
representing the position of the initial state. In the latter 
case we generate the SMA states under < |<S| < 1 by 
parametric transfer of the population from mp = ±1 to 
mp — from the in-phase stationary state represented 
by point (c) in Fig.[4j Our procedure yields: uq > nxot/2 
and n_i = n+i < nxot/4, which is equivalent to moving 
upwards in the diagram Fig. [TTJ This situates the local 
spin of the engineered state inside of the spin sphere, at 
the "equator" , between point (c) and point (b) in Fig. [4] 
Both of the above procedures result in initial states that 
exhibit spin-mixing dynamics along closed phase space 
trajectories. We emphasize that not only the amplitude 
but, most importantly, also the period of the spin-mixing 
oscillations depends on the degree of the population im- 
balance, i.e. the larger the imbalance, the longer the pe- 
riod of the oscillations (cf. Q). 

In order to model the process of controlled manipula- 



10 



tion of the coherent spinor dynamics, we simulate con- 
densate dynamics under the influence of a pulsed mag- 
netic field, as described in Ref. Q. As we show below, 
the resulting management of the relative phase due to 
the quadratic Zeeman effect can lead to the arrest of the 
spin-mixing oscillations in the spatially localized state of 
the spinor condensate described by the SMA. As this ef- 
fect is reversible, it allows flexible manipulations of the 
multi-component gap solitons. 

The application of the magneticthc external magnetic 
field B along the z-axis results in a Zeeman shift of hy- 
pcrfine states: 



t)t 



Aa(|*±l| 2 + |*0| 2 -|* T l| 2 )*±l 



Of 



A Q (|*_i| 2 + |*+i| 2 ) * 



(17) 



where Ej is the Zeeman energy of the jth hypcrfinc 
state in transverse harmonic oscillator energy units Ej = 
Ejtkv±. The energy of an atom in a hypcrfinc state with 
electronic spin 1/2 and the nuclear spin 7 = 3/2 (for 
87 Rb), in an external magnetic field B is given by the 
Breit-Rabi equation [H, 0] : 



E±i 
E 



— T gmB — y/l±a + a 2 

A h f A h f r— — g 
r Vl + a 2 - 



(18) 



Here Ahf is the hyperfine splitting at zero magnetic field 
and a = (gifiiB + gjfi^B)/ Ahf, where «i is the nuclear 
magneton, /ib is the Bohr magneton, gi and <?j are the 
gyro-magnetic ratios of electron and nucleus. The mag- 
nitude of the magnetic field is taken to be B = 1 Gauss. 
For 87 Rb [53, H3: 

A hf = 6834.68261090434(3)MHz, 
gi = 0.9951414(10) x 10" 3 , 
g 3 = 2.00233113(20). 

In order to remove fast phase rotation from the numerics, 
we rescale Eqs. (TU) to Ej(B = 0) = 0. This yields the 
following: E±\ = 4^(1 - Vl ± a + a 2 ) =F giUiB, and 

If the Zeeman energy is much greater than other terms 
in equations (|17[) . the wave functions can be assumed to 
acquire phases Qj(t) = —Ejt/K. The relative phase is 
then given by AO = 28 Q - (6 +1 + 0_ x ) = -AEt/H, where 
AE = 2E — E-i - E + i. From Eqs. JTSJ) it can then 
be seen that due to the quadratic Zeeman effect, the hy- 
perfine states will acquire a change in the relative phase 
AO. For a specific phase change, the duration of the mag- 
netic field pulse ts is obtained as ts = A0/AE, which 
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FIG. 12: Arrest of spin-mixing oscillations shown in Fig. [TT] 
due to quadratic Zeeman effect. The magnetic field of mag- 
nitude B — 1 Gauss is applied over the period of ts = 
(A6)/(AE) ~ 2 (1.4 ms). 
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FIG. 13: Density profiles of the hyperfine components of the 
vector gap soliton corresponding to £fl Ila i = 1000 in Fig ll2l 
Periodic lattice potential (shaded) has been increased by a 
factor of 5 for better visualization. 



is around ts — 2 time units (1.4 ms) for the parameters 
chosen here, and is much smaller than the total duration 
of the evolution. The magnetic field pulse is initiated 
when the population no approaches the value 1/2. The 
resulting halt of the spin mixing oscillations is shown 
in Fig. [12J The soliton maintains its shape through the 
transition to the stationary state, with the final state 
shown in Fig. [T3] 

We point out that the advantage of studying the spin 
mixing dynamics with the gap solitons is twofold. First, 
due to the tight confinement in the lattice, solitons ex- 
hibit very regular and undisturbed oscillations. Secondly, 
the stability of the oscillations, together with the suscep- 
tibility to the magnetic field opens up the possibility for 
controlled spinor gap soliton manipulations. 



IX. CONCLUSIONS 

We have modeled the nonlinear behaviour of spin- 
1 Bosc-Einstcin condensates with ferromagnetic ( 87 Rb) 
and polar ( 23 Na) types of spin-dependent interactions 
loaded into a one-dimensional optical lattice. By per- 
forming numerical analysis of the condensate evolu- 
tion, we established that nonlinear periodic Bloch states 
formed by the condensate at the edge of the Brillouin 
zone are modulationally unstable, irrespective of the type 
of interactions. This instability underpins the existence 
of the multi-component, spatially localized states of the 
condensate inside the band gaps of the linear Bloch- wave 
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spectrum. 

Furthermore, we have demonstrated that the local- 
ization of the spinor matter waves in an optical lattice 
takes the form of vector gap solitons and self-trapped 
waves. Both types of localized states may exist as in- 
phase and out-of-phase stationary spin configurations in- 
dependent of whether the inter-atomic interactions are 
polar or ferromagnetic. The properties of these gap states 
are well described within the single-mode approximation 
frequently employed in the description of spinor systems. 

We have also identified nonlinear stationary states that 
cannot be captured by the single mode approximation. 
These are "hybrid" gap states that are composed of a 
spatially localized soliton-like state in one of the hyper- 
fine components and self-trapped states with a phase 
"kink" in the other components. We have examined 
stability properties of non-SMA states and found that 
off-site states exhibit Joscphson-likc oscillations between 
lattice wells, which eventually lead to a decay of the sta- 
tionary state. Furthermore, we have shown that all the 
non-SMA states exhibit self-magnetisation, an intrinsic 
precession of the local spin. 

Finally, we have shown that localized gap solitons pre- 



pared with a population or phase imbalance, i.e. in a 
non-stationary configuration, may display coherent un- 
damped spin- mixing dynamics. Application of an exter- 
nal magnetic field can lead to the arrest of the inter- 
species population oscillations and rapid transfer to an 
equilibrium state. This opens up possibilities for con- 
trolled manipulation of the spinor condensate as a means 
of quantum state engineering, which is an essential ele- 
ment of quantum computing schemes. 
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